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ABSTRACT 


Finite  amplitude  acoustic  standing  waves  in  a  rec¬ 
tangular  air-filled  cavity  with  various  wedge-shape 
boundary  perturbations  were  studied  both  experimentally 
and  theoretically. 

The  experimental  results  showed  that  geometrical  per¬ 
turbations  alter  the  finite  amplitude  behavior  of  the 
cavity  and  that  the  nature  of  these  changes  are  in  qualita¬ 
tive  agreement  with  the  predictions  of  the  theory. 

However,  quantitative  agreement  was  not  observed  possibly 
because  the  perturbation  chosen  did  not  satisfy  all  the 
assumptions  of  theory. 

There  does  not  seem  to  be  any  fundamental  obstacle 
preventing  the  choice  of  a  perturbation  that  would  allow 
the  theory  to  be  critically  tested. 
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I.  INTRODUCTION 


Coppens  and  Sanders  [1,2]  developed  a  non-linear 
acoustic  model  with  the  dissipative  term  describing  the 
viscous  and  thermal  energy  losses  encountered  at  the 
walls  of  a  rectangular  rigid  cavity. 

Several  researchers  [3,4]  have  examined  the  prob ams 
and  experimental  results  are  in  excellent  agreement  with 
the  theory  except  when  a  degeneracy  exists  in  the  cavity. 

The  purpose  of  this  research  was  to  examine  the 
finite  amplitude  behavior  at  the  cavity  in  the  region 
for  which  perturbation  effects  can  be  accurately  measured 
and  to  compare  the  results  with  the  predictions  of  the 


model. 


II.  BACKGROUND  AND  THEORY 


In  1975  [2]  Coppens  and  Sanders  formulated  a  pertur¬ 
bation  expansion  for  the  non-linear  acoustic  wave  equation 
with  a  dissipative  term  describing  the  thermal  energy  loss 
and  the  non-ideal  nature  of  cavity. 

This  model  predicts  that  when  a  near-degeneracy  exists, 
geometrical  perturbations  provide  a  mechanism  whereby  a 
nearly  degenerate  mode  can  affect  the  finite  amplitude 
behavior  of  the  cavity.  (Most  of  the  following  theoretical 
developments  were  quoted  from  [2].) 

The  non-linear  wave  equation  of  the  viscous  fluid  [2] 


is, 


1=0  D2^’ 


p  C  2 
*o  o 


a 

3t  • 


[<F-> 


y,2  ,  8-1  ,  P 


(1) 


o  _  co'o 

where  X.  is  an  operator  describing  the  physical  processes 
for  absorption  and  dispersion.  The  right  side  term  can 
be  interpreted  as  a  virtual  source  distributed  within  the 
cavity  created  by  the  finite-amplitude  interaction  of  the 
linear  waves. 

The  standing  pressure  wave  in  the  rectangular  cavity 
has  the  form 


p  ~  cos  k  X  cos  ky  Y  cos  k  Z  cos  wt 
X  z 

If  the  cavity  is  being  driven  near  resonance,  the 
non-resonant  terms  contribute  almost  nothing  with  respect 
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to  the  resonant  terms.  So  the  standing  wave  pressure 
has  the  form  of 


P  =  ?  Pn  (2) 

n=l  n 

p 

where  p-  P—  =  MR„  cos  nk  X  cos  nk  Y  cox  n  k  Z  sin  (nwt  +  <j>  ) 
PQCo^  n  x  y  z  yn 

By  substituting  (2)  into  (1) ,  a  set  of  coupled,  non-linear, 
transcendental  are  obtained; 


n 


rCOS  1 

lsin^ 


( 4>  - 

v  yn 


9  )  =  NM8Q 
n  n 


cos  6 


n 


n-1 

lh  S 

j=l 


R. 

3 


n-3 


cos 

sin 


( 1>  •  • ) 

3  n-3 


for  all  n>l,  N  = 


-  £ 
j  =  l 


R  ,  .  R .  {cos}  (<D  .  -  <j> . )  ]  (3) 

n+3  3  sm  Tn+3  3 


where  D  indicates  the  spatial  dimension  (D  =  1,  2,  3). 

If  a  geometrically  perfect  rectangular  cavity  is  driven  at 
frequencies  near  the  resonance  frequency  of  the  (0,1,0)  mode, 
Aqo  =  0  and  only  the "family  members"  (o,n,o)  will  contribute 
to  the  finite  amplitude  behavior. 

In  this  research  we  were  interested  in  the  behavior  when 
a  degeneracy  exists  between  the  020  and  100  modes  with  small 
perturbations  on  the  boundary.  By  using  the  boundary  condi¬ 
tion,  p-n  it  can  be  shown  that 

p  '  p2w  *  EP'  +  E2  P"  +  *  p2w  +  £P'  (“) 


where  P2W  is  the  pressure  of  (0,2,0)  mode,  t  is  the  pertur¬ 
bation  constant  defined  as  the  ratio  of  a  geometric 


t 
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perturbation  magnitude  (A)  to  the  width  (Lx)  and  p  is 
(1,0,0)  mode  component. 

Equation  (4)  can  be  expanded 


P  =  p  _  + 

*o2o 


_  i (2wt  +6  +  i  )  , _ . 

Aoo  S100  P100  00  (5) 


where  S10Q  r  Q1Q0  sin  t1Q0  x 


'100 


(1  +  F 


100 


- — ,  A  is  the 

2  j  oo 


first  coefficient  of  the  Fourier  transform  of 
Equation  (5)  describes  how  much  energy  of  the 
transferred  to  the  (0,2,0)  mode. 


the  perturbation. 
(1,0,0)  can  be 
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III.  DEFINITION  OF  SOME  PARAMETERS 


Several  parameters  were  defined  for  the  purpose  of 
simplifying  the  mathematical  formulation  wi-  lout  loss 
of  physical  contents. 

A.  STRENGTH  PARAMETERS  (STRM) 

STRPM  =  MPQ]_  characterizes  the  strength  of  the 
finite-amplitude  interaction.  Since  the  microphone  sensi¬ 
tivity  (Sm  obtained  with  B&K  4220  pistonphone)  is  known, 

STRPM  s  M8Q1  =  7.07  x  10~3 

JO  \7i  ^  C 

where  M  =  ,  3  a  >s(l  +  p2)  =1.2  with  =£  =  1.4 

P(TSM  cv  v 

M 

P  =  1.293  Kg/m 3 ,  c  =  345  m/sec 

=  RMS  output  voltage  of  the  first  harmonic  component 

B.  FREQUENCY  PARAMETER  (FP) 

2Q  ( f-nf . ) 

Fp  =  _D__ - L- 

nfl 

where  f  is  a  resonant  frequency  of  (0,1,0)  mode,  it  is  the 

Aw 

driving  frequency  it  is  assumed  that  <<1. 

n 

FP  indicates  where  the  driving  frequency  is  with  respect 
to  the  fundamental  resonance  frequency  of  the  cavity.  As 
j FP S  increases,  the  non-linear  effects  become  smaller. 
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C.  HARMONICITY  COEFFICIENT  (E(n) 

r,n1  _  fn-nfl  _  Fpn 
E(n)  -  — —g - ott- 


E(n)  characterizes  how  well  the  modes  of  a  given 
family  are  tuned.  If  E(n)  <4j,  the  corresponding  harmonic 
will  be  strongly  excited. 


IV.  THE  PRESSURE  DISTRIBUTION  WITH  A  WEDGED  PERTURBATION 


y 


The  pressure  distribution  of  (0.2,0)  mode  has  the  form 

p._rt  =  P  cos  y  cos  (2wt  +  d~)  (1) 

020  Ijy  *■ 

The  equation  for  the  perturbed  boundary  is 

X  =  L  {1-  p  [  P  (y-fv  ]  [U  (y — 4^’  “  U(y-Ly)]}  (2) 

x  y  * 

=  L  [ 1+e  f (y ,z)  ] 

X 

Where  e  s  A  ,  f  (y,z)  =-^-(y'|L  )  [U  (y— ^-) -U(y-Ly)  ] . 

bx  y 
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Differentiating  p 


and  f{y,z)  with  respect  to  y,  gives 


020  _2ttP  2  it 

-jy-  =  “jj —  sm  y  cos  <2wt+82)  , 


mi  [U(y._Z)  -  U  (y— Ly )  ] 


■  r~  (y  '  I  V [5(y  ■  “4Z) 


-  6(y-Ly) ] 


Equation  (3)  ;is  substituted  into  the  following  Equation 


3P'  i  _  r  <■  3f  (y,2)  ?p020, 

3x  1  .  xl  3y  3y  Jx*L 


giving 


2  -  J  v 

*  A  sin  =—  y  cos  (2wt  +  ‘U  (y — 

x=Ly  y 


-  u (y-Lv) 


+  (y-|Ly)  n  (y— 4^)  -  '  <y-Ly>  ]  > 


8ttPL 

where  A  =  - ^  . 


Equation  (4)  can  be  expanded  as  a  Pourier  series. 


Hr  lX-L  *  A  COS  (2wt  *  V  Am  cos  F  V 

x  m=0  y 


it'  ■  in  n  i 

+Ji  b"slnSyl 


Only  Aq  needed  to  determine  the  first  order  correction  term 
due  to  the  degerate  between  (0,2,0)  and  (1,0,0)  modes. 

Ao  ’  L  ^  sin  y  (u  '  u  ,y-Ly> 

+  (y  -  |  iy)  [«(y-|ty>  -  5  <y-Ly)  1 >  dy 

_1_ 

2  TT 


So  the  perturbation  correction,  P'  is 


P'  * 


2A  A 

o  _  .  TTX  i(2wt  +  eo 

~  2  Q100  Sin  T100  COS  L  6  2 

L  (— )  x 

x  Co 


':100}  (6 


The  pressure  at  microphone  position,  x=L  ,  can  be  written  as 

X 

following  with  (^) 2  =  (^-)2, 

o  Ly 


P  =  p  +  cn1  a  PFte  f  2wt  +  +  —  £A  P  sin  x  x 

mic  020  ~p  ie  n  tA0  Fioo  s  n  T100  x 


ei(2wt  +  2  + 


(7) 


The  total  pressure  amplitude  at  a  microphone  position  will  be 
I pmic !  =  P/(1+Bcos  T100) 2  +  <B  sinx100)2  (8 

where  B  =  |  e  Aq  Q1q0  sin  t10Q 
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[In  reference  [ 4 ] ,  the  author  made  a  error  in  formulating 
Equation  (2)  ;  he  user’  y  instead  of  (y-aL^)  and  got  an 
incorrect  value  of  A0.] 


V.  APPARATUS 


A.  THE  RECTANGULAR  CAVITY 

The  cavity  is  constructed  from  0.75  inch  aluminum 
(Fig.  1) .  The  interior  of  dimension  of  the  cavity  is 
12.002  inches  long,  2,5002  inches  high  with  a  width  that 
can  be  varied  between  5.50  inches  to  7  inches  in  0.25  inch 
increments.  It  was  sealed  by  a  thin  silicon  grease  layer 
prior  to  assembly.  There  are  three  different  microphone 
positions  for  the  purpose  of  discriminating  the  different 
modes. 

Figure  2  shows  the  pressure  distribution  for  several 
of  the  lower  modes.  According  to  Table  I  (calculated  for 
cavity  dimensions,  12"x6"x2 . 5" )  the  (0,2,0)  and  (1,0,0) 
modes  are  degenerate.  Experimentally  (0,2,0)  mode  was 
about  3  Hz  higher  than  (1,0,0)  mode. 

The  piston  is  set  flush  with  the  bottom  of  the  cavity 
as  near  to  the  wall  Y=0  as  practically  possible,  and  half¬ 
way  between  the  wall  at  X=C  and  X=Lx«  In  this  position 
it  can  excite  the  (0,n,0)  family  of  modes  but  it  cannot 
excite  the  (1,0,0)  mode.  To  drive  the  (1,0,0)  mode  it  is 
necessary  to  insert  an  auxilliary  driver  (ID-30)  at 
Position  A.  A  microphone  at  Position  B  will  sense  the 
pressure  of  all  modes.  It  determines  the  resonance 
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frequencies  and  Q's  of  the  degenerate  modes.  The  microphone 
is  placed  at  Position  A  to  sense  the  (1,0,0)  mode  without 
contamination  from  the  (0,2,0)  and  at  Position  C  to  sense 
the  (0,2,0)  mode  without  contamination  from  the  (1,0,0) 
mode. 

B.  BLOCK  DIAGRAM  (Fig.  4) 

GR  1161-A  coherent  decade  frequency  synthesizer  produces 
a  driving  signal  which  is  precise  to  +0.001  Hz.  This 
signal  enters  a  2120  MB  power  amplifier  which  drives  the 
shaker.  The  piston  motion  was  continuously  monitored  by  an 
Endevco  Model  2215  Accelerometer  mounted  within  the  piston, 
measured  on  a  HP  400D  Voltmeter  and  the  output  of  the 
accelerometer  was  observed  qualitatively  on  a  Model  13D  BR  HP 
oscilloscope.  A  Schlumberger  spectrum  analyzer  was  used  to 
check  for  harmonic  distortion  in  the  piston  motion. 

The  output  of  a  microphone  with  matching  preamplifier 
E&K  2801  fed  into  the  three  devices;  (1)  HP  400D  VTVM  for 
measuring  overall  voltage  leve,  (2)  the  spectrum  analyzer 
for  analyzing  the  wave  component,  and  (3)  two  HP  3 02 A  wave 
analyzers  for  measuring  the  pressure  of  the  first  two 
harmonics. 
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Freauency  (Hz) 


Frequency  (Hz) 


010 

565 

020 

1131 

030 

1697 

100 

1131 

110 

1205 

120 

1599.5 

2039.4 

2529.0 

150 

3032.0 

111 

2980.5 

2757.5 

2926.4 

Cavitv  Orientation,  (b)  nressure  'jodes  for  ( 
For  (6,2,0)  ,  (d)  For  (1,0,0) . 
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VI.  DATA  COLLECTION  PROCEDURE 


Since  the  resonance  frequencies  of  the  cavity  was 
observed  to  drift  with  time,  the  system  was  allowed  warm-up 
for  at  least  one  hour  prior  to  data  collection.  The 
harmonic  content  of  the  accelerometer  output  was  analyzed 
and  the  piston  was  adjusted  until  the  second  harmonic  of 
the  accelerometer  output  was  at  least  50  dB  below  that  of 
the  fundamental  frequency.  Figure  5  shows  typical  results 
for  piston  V2//Vi  re<3^-on  the  finite  amplitude 

measurement.  For  a  driving  voltage  typical  of  those  used 
in  the  experiment,  measurements  were  limited  to  those 
frequencies  for  which  V2/V1<1%. 

There  are  three  measurement  procedures  for  collecting 
data;  (1)  a  pre-run  infinitesimal  amplitude  measurement, 

(2)  the  finite  amplitude  measurement,  and  (3)  a  post-run 
infinitesimal  measurement. 

A.  PRE-RUN  AND  POST-RUN  INFINITESIMAL  MEASUREMENTS 

During  these  measurements,  the  piston  was  driven  so 
that  accelerometer  output  was  less  than  0.1  v.  The  resonant 
frequency  of  the  fundamental  and  the  Q  factor  first  few 
(0,n,0)  in  the  modes  were  determined  from 
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The  resonance  frequency  and  Q  of  (1,0,0)  were  determined 
by  driving  the  cavity  at  port  "A: . 

The  harmonicity  coefficients  were  then  calculated 

from 


E 


n 


Vnfi 


for  E 


ono 


for  E 


100 


The  times  at  which  f  and  fT  were  taken  was  recorded. 

U  Lj 


B.  THE  FINITE  AMPLITUDE  MEASUREMENT 

Throughout  this  measurement  had  to  be  kept  constant 
which  required  adjustment  of  the  driving  voltage  applied 
to  the  piston.  The  HP  302A  wave  analyzer  was  set  to  the 
driving  frequency  and  the  driving  voltage  adjusted  to 
keep  constant  as  the  frequency  was  changed  with  the 
microphone  in  Postion  B.  The  frequency  was  increased 
0.5  Hz  steps  through  +6  Hz  about  the  resonance  frequency 
of  a  fundamental  mode.  At  each  driving  frequency,  the 
harmonic  content  of  the  microphone  onput,  usually  up  to 
the  fourth,  was  measured  with  HP  302A  wave  analyzer  set 
on  AFC  mode  and  the  time  of  that  measurement  recording 
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presented.  The  results  were  presented  as  V^/V^  vs.  the 
frequency  parameter  for  a  given  strength  parameter  and 
perturbation. 

To  determine  the  exact  value  of  the  frequency  parameter, 

the  value  of  f  at  the  instant  of  the  measurement  of  V  was 
r  n 

determined  by  interpolation  of  the  values  of  f  found  in 
the  pre-  and  post-run  procedures.  Figure  6  shows  the 
drift  in  f  is  approximately  linear  in  time. 
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VII.  RESULTS 


Figure  7  shows  the  excellent  agreement  between  theory 
and  experiment  obtained  when  there  are  no  perturbations  so 
that  only  one  family  of  modes  is  excited. 

Some  discrepancies  were  always  observed  when  the  driving 
frequency  was  far  from  the  fundamental  resonance  frequency 
(for  | FP | <2 . )  In  this  region,  the  piston  had  to  be  driven 
hard  to  keep  V,  constant.  This  causes  a  second  harmonic  to 
appear  in  the  piston  waveform  (as  shown  by  Fig.  3) ,  thereby 
introducing  a  linearly-generated  second  harmonic  into  the 
cavity. 

For  the  boundary  perturbation  of  Fig.  3,  the  results 
are  plotted  in  Figs.  8  through  15.  The  thin  line  is  the 
theory  prediction  for  the  perturbation  and  the  thick  line 
includes  the  perturbation  correction.  The  small  circles 
are  experimental  results.  Along  the  edge  of  each  figure 
is  the  information  about  the  cavity  configuration. 

A.  THE  WEDGE  AT  THE  CORNER  OF  THE  CAVITY 

Figures  8  through  11  show  the  results  for  a  wedge  at 
the  corner  of  the  cavity.  The  wedge  used  are  those  of 
Fig.  3  (a) -(d)  respectively. 

B.  THE  WEDGE  AT  T!-7rn  CENTER  OF  THE  LONG  WALL 

Figures  12  through  15  show  the  results  for  the  same 
wedges  (Fig.  3  (a) -(h))  but  now  located  at  the  center  of 


28 


0.177  with 


the  long  wall.  For  Figures  12  and  13,  Aq  = 
DELTA  for  Figure  13  twice  that  for  Figure  12. 
and  15,  Aq  =  0.500  with  DELTA  the  same  as  for 
13. 


For  Figures  14 
Figures  12  and 
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VIII.  CONCLUSIONS 


A.  The  experimental  apparatus  and  procedures  are 
capable  of  providing  data  sufficiently  precise  to  verify 
the  prediction  of  the  theory  in  the  absence  of  geometrical 
perturbation  if  the  range  of  driving  frequencies  is 
restricted  so  that  the  absolute  value  of  the  frequency 
parameter  is  less  than  2. 

B.  To  verify  the  correction  for  a  geometrical  per¬ 
turbation,  it  is  useful  to  have:  (1)  a  sufficiently  large 
correction  to  the  unperturbed  prediction  and  (2)  the  sign 
of  the  correct  must  change  for  an  absolute  value  of  the 
frequency  parameter  less  than  2. 

C.  The  experimental  results  show  that  geometrical 
perturbation  alters  the  finite  amplitude  behavior  of  the 
cavity,  and  that  the  nature  these  change  are  in  qualita¬ 
tive  agreement  with  the  predictions  of  the  theory.  However, 
quantitative  agreement  was  not  observed. 

D.  Sources  of  the  difficulties  in  obtaining  good 
agreement  might  be;  (1)  the  inability  to  experimentally 
satisfy  conditions  of  B  above,  (2)  the  DELTA  function 
arising  from  the  perturbations  used  may  cause  trouble  in 
the  theory  (Also,  they  may  be  difficult  to  accurately 
reproduce  experimentally.),  and  (3)  the  higher  order 
terms  neglected  in  developing  the  theory  may  not  all  be 


small.  For  example,  terms  of  order  higher  than  first 
order  in  Equation  (4)  of  Section  II  and  terms  of  order 
higher  than  Aq  and  all  bm  in  Equation  (5)  of  Secv.ion  IV. 

E.  The  theory  for  the  perturbed  cavity  can  probably 
be  critically  tested,  if  a  perturbation  shape  is  used 
that  avoids  the  difficulties  of  D(2)  and  satisfies  the 
condition  of  B.  At  this  time  it  is  not  obvious  that 
such  a  perturbation  can  be  found,  but  there  seems  to  be 
no  fundamental  obstacle  to  find  a  suitable  perturbation. 
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TABLE  I  (Fig.  7) 


mode 

time 

f 

u 

fL 

f 

r 

Q 

E 

010 

0 

565.77 

562.88 

564.33 

195.40 

o  _2 

020 

3 

1131.84 

1127.79 

1129.81 

278.69 

0.10x10  z 

030 

6 

1698.27 

1693.52 

1695.90 

356.51 

0.17xl0'2 

100 

9 

1131.87 

1127.59 

1129.73 

264.20 

-0.74xl0“4 

ii 

lldB 

time  f 

V2 

V3  V4 

P2/Pl 

p3/Pl  1 

12 

558.5 

-50.0 

-70.5 

-70.5 

13 

559.0 

-47.9 

-72.0 

-69.0 

14 

559.5 

-46.9 

-73.8 

-67.8 

15 

560.0 

-46.0 

-72.6 

-67.5 

16 

560.5 

-44.8 

-71.2 

-69.7 

17 

561.0 

-43.7 

-69.7 

-73.5 

18 

561.5 

-42.4 

-67.8 

-76.5 

19 

562.0 

-41.0 

-65.7 

-75.9 

20 

562.5 

-39.4 

-63.2 

-82.5 

21 

565.0 

-37.8 

-60.3 

-85.0 

22 

563.5 

-35.9 

-56.6 

-74.9 

23 

564.0 

-34.0 

-52.2 

-67.5 

24 

564.5 

-32.2 

-47.6 

-60.6 

25 

565.0 

-31.1 

-45.0 

-54.5 

26 

565.5 

-31.6 

-43.9 

-52.9 

27 

566.0 

-33.3 

-46.7 

-56.0 

28 

566.5 

-35.3 

-51.0 

-63.2 

29 

567.0 

-37.6 

-55.5 

-70.8 

30 

567.5 

-39.8 

-60.0 

-73.6 

31 

568.0 

-42.1 

-64.3 

-74.7 

32 

568.5 

-44.1 

-68.5 

-74.2 

33 

569.0 

-46.2 

-69.9 

34 

569.5 

-47.2 

-68.7 

-67.0 

35 

570.0 

-48.0 

-65.4 

-67.0 

36 

570.5 

-48.9 

-65.1 

-67.2 

37 

571.0 

-49.8 

-65.2 

-67.5 

38 

571.5 

-50.8 

-65.4 

-68.2 

39 

572.0 

-51.7 

-65.8 

-68.6 

40 

572.5 

-52.5 

-66.3 

-67.4 

41 

573.0 

-54.7 

-65.1 

-66.7 

0.0113 

0.0011 

0.0143 

0.0009 

0.0160 

0.0119 

0.0204 

0.0011 

0.0233 

0.0012 

0.0269 

0.0014 

0.0316 

0.0018 

0.0380 

0.0025 

0.0457 

0.0034 

0.0569 

0.0052 

0.0708 

0.0007 

0.0876 

0.0148 

0.0994 

0.0200 

0. 0933 

0.0226 

0.0772 

0.0164 

0.0610 

-0.0100 

0.0484 

0.0060 

0.0365 

0.0036 

0.0279 

0.0022 

0.0221 

0.0013 

0.0174 

0.0011 

0.0156 

0.0013 

0.0141 

0.0019 

0.0127 

0.0020 

0.0115 

0.0019 

0.0103 

0.0019 

0.0093 

0.0018 

0.0084 

0.0017 

0.0065 

0.0020 

0.0011 

-4.19 

0.0013 

-3.85 

0.0014 

-3.50 

0.0015 

-3.15 

0.0012 

-2.81 

-2.47 

-2.12 

-1.78 

-1.43 

-1.09 

— 

-0.74 

0.0015 

-0.40 

0.0033 

-0.05 

0.0067 

0.29 

0.0081 

0.64 

0.0056 

0.98 

0.0025 

1.33 

0.0010 

i.67 

2.02 

2.36 

2.71 

3.05 

0.0016 

3.40 

0.0010 

3.14 

0.0016 

4.09 

0.0015 

4.43 

0.0014 

4.78 

0.0013 

5.12 

0.0015 

3.47 

0.0016 

5.81 

mode 

time 

f 

u 

fL 

f 

r 

Q 

E 

010 

44 

566.28 

563.37 

564.82 

194.03 

0 

020 

47 

1132.74 

1128.67 

1130.71 

277.82 

0.94xl0”3 

030 

50 

1699.77 

1694.97 

1697.37 

353.84 

0. 17x1 0-3 

100 

53 

1132.80 

1128.48 

1130.64 

262.03 

-0.61X10-4 
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mode  time 


E 


010 

0 

568.27 

565.28 

566. 

77 

189.81 

0 

020 

3 

1132.37 

1127.62 

1129. 

99 

238.25 

-0.31x10- 

-2 

030 

6 

1701.89 

1696.79 

1699. 

34 

333.14 

-0.57x10* 

-3 

100 

9 

1132.44 

1127.71 

1130. 

07 

238.76 

0.69x10' 

-4 

V^-lldB 

time 

f 

V2 

V3 

V4 

V*1 

P3^1 

V*1 

FP 

12 

561.3 

-44.3 

-69.9 

-71.3 

0.0216 

0.0011 

-3.95 

13 

561.8 

-43.2 

-69.2 

-72.2 

0.0247 

0.0012 

-3.61 

14 

562.3 

-41.7 

-68.3 

-74.0 

0.0293 

0.0014 

-3.26 

15 

562.8 

-40.2 

-66.8 

-80.0 

0.0347 

0.0016 

-2.92 

16 

563.3 

-38.5 

-65.2 

0. 0422 

0.0020 

-2.58 

17 

563.8 

-36.4 

-62.9 

0.0537 

0.0025 

-2.24 

18 

564.3 

-34.3 

-60.3 

0.0684 

0.0034 

-1.90 

19 

564.8 

-32.8 

-57.2 

0.0813 

0.0049 

-1.56 

20 

565.3 

-32.5 

-55.0 

-72.0 

0.0841 

0.0063 

0.0009 

-1.22 

21 

565.8 

-33.5 

-53.5 

-67.8 

0.0750 

0.0075 

0.0014 

-0.87 

22 

566.3 

-35.1 

-51.8 

-62.5 

0. 0624 

0.0092 

0.0024 

-0.53 

23 

566.8 

-36.9 

-50.9 

-60.9 

0.0507 

0.1002 

0.0032 

-0.19 

24 

567.3 

-37.9 

-51.2 

-60.8 

0.0452 

0.0098 

0.0032 

0.15 

25 

567.8 

-38.6 

-52.8 

-64.6 

0.0419 

0,0082 

0.0021 

0.49 

26 

568.3 

-39.2 

-55.3 

-69.2 

0.0391 

0.0061 

0.0072 

0.83 

27 

568.8 

-40.3 

-58.7 

-75.0 

0.0343 

0.0041 

1.78 

28 

569.3 

-41.5 

-61.7 

-78.5 

0.0300 

0.0029 

1.52 

29 

569.8 

-42.4 

-64.3 

0.0271 

0.0022 

1.86 

30 

570.3 

-44.3 

-67.3 

0.0218 

0.0015 

2.20 

31 

570.8 

-45.5 

-70.8 

0.0188 

0.0010 

2.54 

32 

571.3 

-46.9 

-75.0 

0.0160 

0. 0006 

2.88 

33 

571.8 

-48.1 

-76.9 

0.0140 

0.0005 

3.22 

34 

572.3 

-49.3 

-75.2 

-69.1 

0.0122 

0.0006  , 

0.0012 

3.57 

35 

572.8 

-50.6 

-73.5 

-68.0 

0.0105 

0.0007 

0.0014 

3.91 

36 

573.3 

-51.0 

-71.7 

-66.5 

0.0100 

0.0009 

0.0017 

4.'’' 

37 

573.8 

-51.8 

-71.1 

-66.0 

0.0092 

0.0010 

0.0018 

4.o. 

38 

574.3 

-52.5 

-70.3 

-65.7 

0.0084 

0.0011 

0.0019 

4.93 

39 

574.8 

053.8 

-72.3 

-66.7 

0.0073 

0.009 

0.0016 

5.27 

mode 

time 

f 

u 

f2 

f 

r 

Q 

E 

010 

42 

568.82 

565.95 

567.38 

197.49 

0 

020 

45 

1133.50 

1128.70 

1131.10 

235.40 

-0. 32x1 O-2 

030 

48 

1703.98 

1698.69 

1701.34 

321.25 

-0.48xl0"3 

100 

51 

1133.85 

1129.10 

1131.47 

238.00 

0.33xl0"3 
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TABLE  III  (Fig.  9) 


mode 

time 

f 

u 

f2 

f 

r 

Q 

E 

010 

0 

567.85 

565.01 

566.43 

199.52 

0 

020 

3 

1128.61 

1124.81 

1126.71 

296.58 

-0.54xl0-2 

030 

6 

1697.95 

1693.16 

1695.39 

354.31 

-0. 23xl0~2 

100 

9 

1128.77 

1124.87 

1126. 82 

288.63 

-0.96xl0-4 

V^=-lldB  time 

f  V2 

V3 

v4  p2/P1  p3/p1 

VP1  ] 

12 

561.0 

-37.6 

13 

561.5 

-36.2 

14 

562.0 

-34.3 

15 

562.5 

-32.2 

16 

563.0 

-30.2 

17 

563.5 

-29.5 

18 

564.0 

-30.8 

19 

564.5 

-33.3 

20 

565.0 

-35.5 

21 

565.5 

-37.3 

22 

566.0 

-38.9 

23 

566.5 

-40.3 

24 

567.0 

-41.7 

25 

567.5 

-42.5 

26 

568.0 

-42.6 

27 

568.5 

-43.1 

28 

569.0 

-43.3 

29 

569.5 

-43.8 

30 

570.0 

-44.5 

31 

570.5 

-45.2 

32 

571.0 

-45.6 

33 

571.5 

-46.2 

34 

572.0 

-46.8 

35 

572.5 

-47.5 

36 

573.0 

-47.5 

37 

573.5 

-48.2 

-66.2  0.0468 

-64.8  0.0550 

-62.8  0.0684 

-60.3  0.0876 

-57.5  -76.0  0.1096 

-55.6  -73.4  0.1189 

-55.2  -74.0  0.1023 

-55.5  -74.7  0.0772 

-55.8  -71.8  0.0599 

-58.0  -70.2  0.0484 

-60.5  -71.4  0.0403 

-61.5  -76.0  0.0345 

-61.7  -80.0  0.0293 

-61.3  -81.3  0.0266 

-60.8  -80.3  0.0263 

-62.2  0.0250 

-63.7  0.0243 

-65.9  0.0230 

-68.0  0.0213 

-69.6  0.0196 

-70.9 - 0.0186 

-71.7  0.0174 

-72.7  0.0162 

-73.3  0.0150 

-73.5  0.0150 

-77.0 - 0.0139 


0.0017  -3.84 

0.0021  -  -3.49 

0.0026  - -3.14 

0.0034  - -2.79 

0.0047  0.0006  -2.44 

0.0059  0.0008  -2.09 

0.0062  0.0007  -1.74 

0.0060  0.0007  -1.39 

0.0058  0.0009  -1.04 

0.0045  0.0011  -0.67 

0.0033  0.0010  -0.34 

0.0030  0.0006  0.00 


0.0029  0.36 

0.0031  -  0.71 

0.0032  1.06 

0.0028  -  1.41 

0.0023  -  1.76 

0.0018  -  2.11 

0.0014 - 2.46 

0.0012 - 2.81 

0.0010 - 3.16 

0.0009  -  3.51 

- - 3.86 

- 4.21 

- 4.56 

- 4.91 


mode 

time 

f 

u 

fL 

010 

40 

567.99 

565.11 

020 

43 

1128.94 

1125.12 

030 

46 

1698.43 

1693.67 

100 

49 

1129.10 

1125.21 

f 

r 

Q 

E 

566.55 

196.79 

0 

-2 

1127.03 

295.42 

-0.54x10 

-2 

1696.05 

355.76 

-0.21x10 

-3 

1127.15 

289.76 

0.11x10 

45 


010 

0 

566.19 

563 

.20 

564.70 

189.30 

0 

020 

3 

1128.86 

1124 

.95 

1126.91 

287.70 

-0.2x10' 

-2 

030 

6 

1696.18 

1691 

.38 

1693.78 

353.24 

-0. 18xl0~3 

100 

9 

1128.46 

1124 

.62 

1126.54 

293.29 

-0.32xl0~3 

V^-lldB 

time 

f 

V2 

V3 

V4 

^2^1 

V*i 

V*1 

FP 

12 

558.7 

-42.8 

-69.8 

0.0257 

0.0011 

-4.04 

13 

559.2 

-41.8 

-68.7 

0.0288 

0.0013 

-  ■?  ->l 

14 

559.7 

-40.8 

-67.9 

0.0324 

0.0014 

*o.  37 

15 

560.2 

-39.7 

-66.9 

0.0367 

0.0016 

-2  04 

16 

560.7 

-38.3 

-65.5 

0. 0434 

0.0019 

-2.71 

17 

561.2 

-36.7 

-63.5 

— 

0.0519 

0.0024 

-2.38 

18 

561.7 

-34.9 

-61.0 

0.0638 

0.0032 

-2.04 

19 

562.2 

-33.0 

-57.8 

-76.3 

0.9794 

0.0046 

-1.71 

20 

562.7 

-30.7 

-53.7 

-70.7 

0.1041 

0.0074 

0.0010 

-1.38 

21 

563.2 

-28.9 

-50.2 

-65.3 

0.1274 

0.0110 

0.0019 

-1.04 

22 

563.7 

-28.5 

-47.5 

-61.4 

0.1334 

0.0150 

0.0030 

-0.71 

23 

564.2 

-30.4 

047.2 

-59.0 

0.1078 

0.0155 

0.0040 

-0.38 

24 

564.7 

-32.2 

-48.5 

-58.5 

0.0873 

0.0133 

0.0042 

-0.05 

25 

565.2 

-33.7 

-52.2 

-62.7 

0.0733 

0.0087 

0.0026 

0.29 

26 

565.7 

-35.6 

-56.4 

-70.07 

'  0.0592 

0.0054 

0.0010 

0.62 

27 

566.2 

-37.0 

-58.4 

-74.7 

0.0501 

0.0043 

0.95 

28 

566.7 

-39.6 

-59.7 

-77.2 

0.0314 

0.0037 

1.29 

29 

567.2 

-42.1 

-61.5 

0.0279 

0.0030 

1.62 

30 

567.7 

-43.7 

-63.7 

0.0232 

0.0023 

1.95 

31 

568.2 

-44.5 

-65.5 

0.0221 

0.0019 

— 

2.29 

32 

568.7 

-45.1 

-66.9 

0.0197 

0.0016 

2.62 

33 

569.2 

-45.6 

-68.0 

0.0186 

0.0014 

2.95 

34 

569.7 

-46.1 

-69.0 

0.0176 

0.0013 

3.29 

35 

570.2 

-46.6 

-69.8 

— 

0.0167 

0.0012 

3.62 

36 

570.7 

-47.0 

-70.4 

0.0158 

0.0011 

3.95 

37 

571.2 

-47.2 

-70.6 

0.0155 

0.0011 

4.28 

38 

571.7 

-47.7 

-71.3 

— 

0.0146 

0.0010 

— 

4.62 

39 

572.2 

-48.6 

-63.9 

0.0132 

0.0023 

4.95 

mode 

time 

f 

fT 

f 

Q 

E 

___ 

u 

L 

r 

010 

42 

566.35 

563.33 

564.84 

186.79 

0 

020 

45 

1129.43 

1125.52 

1127.48 

1287.84 

-0. 19xl0~2 

030 

48 

1697.04 

1692.26 

1694.65 

353.94 

-0. 7 8x1 0-4 

100 

51 

1129.12 

1125.28 

1127.20 

294.16 

-0.24xl0-3 

46 


TABLE  V  (Fig.  11) 


mode 

time 

fu 

f  L 

010 

0 

566.90 

563.86 

020 

3 

1127.19 

1123.28 

030 

6 

1695.31 

1690.47 

100 

9 

1126.76 

1122.90 

V1=-lldB 

time 

f  V2 

V3 

f 

r 

Q 

E 

565.38 

185.80 

0 

-2 

1125.23 

288.08 

-0.49x10 

1692.89 

349.27 

-0.19xlQ-2 

-3 

1124.83 

291.41 

-0.36x10 

V4  P2/P1  P3/P!  P4/P1  FP 


12 

559.4 

-39.5 

-66.9 

13 

559.9 

-37.9 

-65.3 

14 

560.4 

-36.4 

-63.9 

15 

560.9 

-34.5 

-61.8 

16 

561.4 

-32.4 

-59.1 

17 

561.9 

-30.3 

-55.7 

18 

562.4 

-28.7 

-55.6 

19 

562.9 

-28.9 

-50.8 

20 
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4.  Professor  James  V.  Sanders,  Code  61Sd 
Department  of  Physics  and  Chemistry 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Professor  Alan  B.  Coppens ,  Code  61Cz 
Department  of  Physics  and  Chemistry 
Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Dr.  D.  T.  Blackstock 
Applied  Research  Laboratory 
University  of  Texas 
Austin,  Texas  78712 

7.  Library,  Naval  Academy 
Chin-Hae,  Kyung  Sang  Nam  Do 
KOREA 

8.  Dr.  Lee,  C.H. 

Add,  Chin-Hae,  Kyung  Sang  Nam  Do 
KOREA 

9.  LT  Ender  Kuntgal 
Akatlar,  Zeytinoglu  Caddesl 
Satak  Ap.  240/16 

Etiler,  Istanbul,  TURKEY 

10.  LT  Mehmet  Aydin 
Zincirlikuyu  Caddesl  No-99 
Kasimpasa,  Istanbul,  TURKEY 

11  .  LTC  Ilbok  Joung 

Anti-Submarine  Bureau 
Navy  Headquarters,  Seoul 
KOREA 
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